Pairing of a harmonically trapped fermionic Tonks- Girardeau gas 
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The fermiorric Tonks-Girardeau (FTG) gas is a one-dimensional spin-polarized Fermi gas with in- 
finitely strong attractive zero-range odd- wave interactions, arising from a confinement-induced reso- 
nance reachable via a three-dimensional p-wave Feshbach resonance. We investigate the off-diagonal 
long-range order (ODLRO) of the FTG gas subjected to a longitudinal harmonic confinement by 
analyzing the two-particle reduced density matrix for which we derive a closed-form expression. 
Using a variational approach and numerical diagonalization we find that the largest eigenvalue of 
the two-body density matrix is of order N/2, where A'^ is the total particle number, and hence a 
partial ODLRO is present for a FTG gas in the trap. 

PACS numbers; 03.75.-b,05.30.Jp 



I. INTRODUCTION 

Low-dimensional systems display unusual and strik- 
ing features as compared to their three-dimensional (3D) 
counterparts, among which are the enhanced effects of 
the interactions and the presence of large fluctuations 
which are responsible for the failure of mean-field ap- 
proaches. In addition, especially in the one-dimensional 
(ID) case, it is possible to find exact solutions which 
greatly help the progress in our understanding of com- 
plex, strongly-interacting many-body systems. 

Low-dimensional atomic quantum gases are one of 
the frontiers of the current theoretical and experi- 
mental investigations. Bosonic and fermionic atomic 
gases constrained to a quasi- ID geometry have al- 
ready been realized experimentally by trapping atoms 
in two-dimensional optical lattices P, In the 

case of bosons the strongly-repulsive (known as Tonks- 
Girardeau) regime has been achieved 0, H . This regime 
has been fully understood thanks to the knowledge of the 
exact many-body wavefunction through the mapping of 
the strongly repulsive, impenetrable bosons onto an ideal 
as of fermions subjected to the same external potential 
Quite remarkably, the mapping does not hold only 
for the ground state, but may be extended to treat time- 
dependent phenomena and systems at finite temperature 

In this paper we focus on the model of spin-polarized 
fermions interacting via strongly attractive odd- wave in- 
teractions, which are the ID analogue of p- wave interac- 
tions in 3D. This regime might be experimentally reach- 
able by exploiting the so-called confinement-induced res- 
onances (CIR), which permit to tune the ID coupling 
constant via a 3D Feshbach resonance 0. In the limit 
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of iirfinitely strong attractions, known as the fermionic 
Tonks-Girardeau (FTG) regime, the many-body wave- 
function is known exactly through an inverse Fermi-Bose 
mapping which allows to express the fermionic wavefunc- 
tion in terms of the one of an ideal Bose gas 

UME^ At 

the resonance point we can thus quantitatively address 
the question of what the structure of the ground state 
is, and in particular whether the fermions are paired. 
This contributes to the understanding of the intermedi- 
ate, strongly-interacting region in the ID equivalent of 
the BCS-BEC crossover for p-wave fermions. 

In a previous work |l2l | we have studied the issue of 
pairing for a homogeneous FTG gas in the thermody- 
namic limit. We have found that there is off-diagonal 
long-range order (ODLRO) in the reduced two-body den- 
sity matrix, indicating a paired state. We consider here 
the experimentally-relevant case of a gas with a finite 
number of particles and subjected to a longitudinal har- 
monic confinement. We derive an an exact analytic ex- 
pression for the two-body density matrix of the trapped 
gas and both by variational estimates and by numeri- 
cal diagonalization we show that also in the presence of 
the trap ODLRO persists, yielding a complex, partially 
paired quantum state. 



II. FERMIONIC TONKS-GIRARDEAU GAS IN 
A LONGITUDINAL HARMONIC TRAP 



We consider an atomic spin-polarized Fermi gas con- 
fined by a tight atom waveguide of which it occupies the 
transverse ground state. We shall henceforth assume that 
the level spacing in the transverse direction is much larger 
than all the relevant energy scales in the problem (such as 
chemical potential, temperature, etc.), so that the gas is 
effectively one-dimensional. We also consider the case of 
a gas subjected to a much weaker longitudinal harmonic 
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confinement of frequency uj. The ID Hamiltonian is 
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(1) 

where vf^^^ is a short-range attractive two-body interac- 
tion which is specified here below. Since the spatial wave 
function is antisymmetric due to the spin polarization, 
there is no zero-range s-waye (delta function) interaction, 
but it has been shown j^, ^3 Ol that a strong, attrac- 
tive, and short-range odd-wave interaction (a ID analog 
of 3D p-wave interactions) occurs in the neighborhood of 
the CIR. Such an interaction can be expressed through a 
contact condition for the relative wavefunction iPf{x) of 
each pair of fermions as [l^ 



4>f{x) dx 



2h^ 



(2) 



where gfjj is the ID fermionic coupling constant, which 
can be expressed in terms of the 3D p-wave scattering 
volume 0. The condition ^ above, together with the 
antisymmetry condition ipF^x < 0) = —ipFix > 0) on 
the relative wavefunction leads to a solution for ipFix) 
which is discontinuous at the contact point x = but 
has a continuous first derivative. 

The FTG gas corresponds to the negative side of the 
CIR, that is, the case where gf^ — oo; in the FTG 
limit the first derivative of the relative wavefunction 
tpF{x) equals to zero at a; — 0. For the sake of il- 
lustration, let us consider first two particles under har- 
monic confinement. In the presence of such a confine- 
ment the Schroedinger equation can be solved exactly 
for any value of the coupling strength |Q| . The rela- 
tive wavefunction ipFix) in the FTG limit takes the value 

i^F{x > 0) = (27r)-i/4a:o//2e-'3'/'* and i^pix < 0) = 
~{2Tr)^^/'^Xo}J^e^'^^^'^, where we have set Q — x/xosc 
and Xosc = \/ h/muj. 

The solution for N — 2 can be extended to arbitrary 
particle numbers N, so that the exact fermionic TG gas 
ground-state wavefunction is 0, ^| 

N 

'^pixi,- ■ ■ ,xn) = A{xi, ■ ■ ■ ,xn) Y\_ ^o{xj) (3) 



with where A{xi,--- ,xn) = Ili<j<e<N '^S^i^i ~ ^j) 
is the "unit antisymmetric function" employed in the 
original discovery of fermionization 5] and (t>o{x) = 
7r~^/*XosV^e~'5 /2 ig ^]^g ground-state orbital of the lon- 
gitudinal harmonic confinement. Hence, the FTG gas is 
mapped through the function A to the ground state of 
an ideal Bose gas under harmonic confinement, of which 
it shares all the properties that do not depend on the 
sign of the many-body wavefunction, such as the density 
profile and the spectrum of collective excitations. 



III. OFF-DIAGONAL LONG-RANGE ORDER 
FOR A TRAPPED FTG GAS 

In order to explore the pairing properties of the FTG 
gas we study the two-body reduced density matrix, de- 
fined as 

P2{xi,X2;x[,x'2) = N{N-1) / •^f{xi,X2--- ,xn) 



X •^}.{x[,X2,X3,- ■ ■ ,XN)dx3 ■ ■ ■ dxN- 

(4) 

As it was discussed by Yang ^ , the criterion for "super- 
conductive" off-diagonal long-range order (ODLRO) in a 
trapped system is that the largest eigenvalue Ai of the 
two-body density matrix is of the order of the number 
of particles N, i.e. Ai = aN, with < a < 1 being the 
pair-condensate fraction. 

Using the exact form (jS)) of the many-body wavefunc- 
tion, the integration over N — 2 variables in Eq. can 
be explicitly performed and we obtain an analytic expres- 
sion for the two-body density matrix, which reads 

P2{xi,X2;x[,X2) = N{N - l)sgn(a::i - X2)(t)o{xi)(l)o{x2) 
X sgii{x[~x'2)Mx'i)Mx'2)[GixuX2;x[,X2)]'^-^, (5) 

where G{xi,X2;x[,x'2) = 1 + crf(?;i) - erf(y2) + erf(y3) - 
erf(y4), and yi < y2 < Vs < Vi are (Qi, Q2, Qi, Q2) in 
ascending order, and Qi = Xi/xosc- 

The eigenvalues Xj and eigenfunctions uj of the two- 
body density matrix are the solutions of the integral 
equation 

dx'i dx'2P2ixi, X2', x'l, x'2)Uj{x[, x'2) = XjUj{xi,X2)- 

(6) 

If = 2 then G{xi,X2',x'i,X2) = 1, the two-body den- 
sity matrix separates trivially and has only one nonzero 
eigenvalue Ai = iV = 2, with eigenfunction ui(xi,X2) — 
sgn(xi — a;2)^o(a^i)^o(a^2)- For = 3 we have found 
by numerical integration that the ansatz ui{xi,X2) = 
C sgn(a;i - a;2)<^o(a;i)0o(a;2)[l - |erf(Qi) - erf{Q2)\] sat- 
isfies the above eigenvalue equation to within 0.1% with 
eigenvalue Ai = 2 (the same as for N = 2), and we con- 
jecture that these expressions are exact. For iV > 4 we 
have numerically solved the eigenvalue equation © by di- 
agonalization. The results for the largest eigenvalue are 
displayed in Fig. ^ for = 4 to 8. The obtained values 
agree with the analytical expression Ai = {N + l)/2 to 
within the computational accuracy of a few percent |l6j | , 
indicating that in the trap a partial ODLRO is present. 
The eigenfunctions corresponding to the largest eigenval- 
ues obtained by numerical diagonalization are shown in 
Fig. 121 along with the analytical eigenfunctions for = 2 
and 3. 

We have also calculated lower bounds to Ai for A^ = 4 
to 8 from the inequality 



Ai > 



/ Utriali'X.)p2(K,X')utriali'X.')dXdX.' 

Jula^)dX 



(7) 
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FIG. 1: Largest eigenvalue Ai of the two-body density matrix 
for a harmonically trapped fermionic TG gas at various par- 
ticle numbers A*' from numerical diagonalization (stars) and 
the variational lower-bound estimate (diamonds). 




FIG. 2: Eigenfunctions of the two-body density matrix of 
a harmonically trapped fermionic TG gas belonging to the 
largest eigenvalue Ai, for A'' = 2 to 7 from top to bottom 
(when x\ > 0). 



where X is a shorthand notation for (a;i, X2)- In particu- 
lar, we have chosen the ansatz utriai{xi, X2) = C sgn(a;i — 
X2)Mxi)Mx2)[l - |erf(gi) - erf(Q2)|]^-'. As in is il- 
lustrated in Fig. n the results are very close to the nu- 
merical eigenvalues, indicating that the ansatz is a good 
guess for the exact solution. There is however some small 
difference, as is shown in Fig. 13 by comparing the above 
ansatz with the numerically determined eigenfunctions in 
the case iV = 4. 

Finally, it is interesting to estimate the behaviour 
of the two-body density matrix in the thermodynamic 
limit, defined as the limit iV — > 00, w — s- keeping 
the central density no — ^'t%{^) = N ^/mLu/nh con- 
stant. In that case Eq. Q takes the same form as in 
the homogeneous system [IjI and hence it follows that 
P2{xi,X2',x'i,x'2) = \iUi{xi,X2)ui{x'i,x'2) (apart from 
terms decaying exponentially with pair separation) , with 
Ai = N/2 and ui{xi,X2) oc sgn(xi - a;2)e-2"»i^i-^^l . 
Hence, the estimate in the thermodynamic limit agrees 




FIG. 3: Comparison of maximal eigenfunction m = mi of the 
two-body density matrix (upper curve when x\ > 0) and vari- 
ational ansatz u = utriai (lower curve when xi > 0) for the 
case A'^ = 4. In the case A'' = 3 (not plotted) the two curves 
(on the same scale) would be indistinguishable. 



to order 1/iV with the numerical solution. 



IV. SUMMARY AND CONCLUDING 
REMARKS 



In summary, by estimating the largest eigenvalue of 
the two-body reduced density matrix we have investi- 
gated the possibility of off-diagonal long-range order for 
a fermionic Tonks- Girardeau gas subjected to a longi- 
tudinal harmonic confinement. The result of variational 
approach and numerical diagonalization together with an 
estimate in the thermodynamic limit is that the largest 
eigenvalue Ai is of order N/2, where N is the parti- 
cle number, and hence we find a partial ODLRO. The 
value Ai ~ N/2 might also be interpreted as the fraction 
of pairs which are Bose-Einstein condensed in the one- 
dimensional equivalent of the BCS to EEC crossover for 
p-wave fermions. The emerging picture is the one of a 
partially paired quantum state, where quantum fluctua- 
tions play a major role in depleting the "EEC" of pairs. 
The consequences of pairing remain to be explored. 
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